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We present the self-consistent Bogoliubov-de Gennes calculations of an stt pairing state of two 
band superconductivity as a model for the FeAs superconductors. The sn state is an s-wave pairing 
state with an internal tt phase, that is, nodeless gaps on each band but with the opposite sign. 
The novel features of this state are investigated by calculating the local density of states of the tt 
phase superconductor/normal metal bilayers. Because of the sign reversal between the two conden- 
sates, the zero bias conductance peak appears as observed in tunneling spectroscopy experiments on 
FeAs superconductors. This eliminates the major obstacle to establish the sit state as the pairing 
symmetry of the FeAs superconductors. 

PACS numbers: PACS: 74.20.Rp, 74.20.Mn, 74.45. +c, 74.50.-l-r 



Introduction - The recent discovery of the iron based 
pnictides superconductors generated enormous interests 
in the community [l|, 0, S I3j @| ■ Like the cuprate super- 
conductors, the pnictides are highly two-dimensional, the 
undoped parent materials are in the antifcrromagnetic 
state, and superconductivity emerges when the antiferro- 
magnetism is suppressed |(|] upon doping either electrons 
or holes into the FeAs planes. The urgent question 
is the pairing mechanism of the FeAs superconductors. 
This will hopefully pave the way to coherently under- 
stand the still elusive pairing mechanism of the high tem- 
perature superconductivity. The first step towards this 
goal is to establish the pairing symmetry of the FeAs 
superconductors. 

Conflicting experimental results are being reported 
with regard to the orbital pairing symmetry. Some ex- 
periments are unambiguously seeing a fully gapped node- 
less superconducting state, while others point towards a 
gap with nodes: The several independent ARPES experi- 
ments on single crystals reported the full gaps around the 
hole Fermi surface although the gap feature is less clear 
around the electron Fermi surface 

i, S, US]- The tem- 
perature dependence of the penetration depth strongly 
suggests a nodeless gap structure [llj . The infrared spec- 
troscopy observed the superconductivity induced features 
which were best described in terms of s-wave pairing [l^ . 

The evidences for nodes are also accumulating: The 
nuclear spin-lattice relaxation rate l/Ti clearly showed 
no coherence peak and ~ dependence below 
U) 15, Also, the resonant spin excitations were 



glet pairing and multiple gaps seems convincing from 
the Knight shift and other thermodynamic measurements 



observed by neutron scattering experiments with the mo- 
mentum transfer corresponding to the magnetic order 
ajid the energy transfer consistent with the Tc scaling 
\it\ . These observations strongly suggest a gap with 
nodes. On the other hand, the evidence for the sin- 



The various tunneling spectroscopy, however, seems 
confusing. The Andreev spectroscopy reoprted by Chen 
and others [3] was well fitted by the s-wave Blonder- 
Tinkham-Klapwijk formalism [l9|, although they ob- 
served zero bias peak (ZBP) which seemed difficult to 
understand within the s-wave pairing formulation. Other 
point contact and tunneling spectroscopy measurements 
also reported the ZBP which strongly suggests a sign 
change of the pairing order parameter 

Many ideas have been put forward to understand the 
seemingly confiicting experimental observations on the 
FeAs materials. Among them, particularly appealing is 
the sign changing s-wave pairing state as advanced by 
Mazin and coworkers [IS]- It seems to be able to explain 
the experimental observations indicating the full gap as 
well as a gap with nodes [25|, was noticed early 

on that there exists this typ e of solution to a multi-band 
BCS gap equation [28, |29|. It is quite exciting that it 
seems to be actually realized in the pnictides. We adapt 
this proposal in this Letter, and call it "stt" pairing state 
because of its natural connection with other tt states. It 
is interesting to note the similarity between the inter- 
nal TT state of the stt pairing and the "external" tt state 
of the superconductor/ferromagnet (S/F) bilayers. It is 
well established that a Cooper pair in the F side of an 
S/F bilayer picks up a non-zero center of mass momen- 
tum and the singlet pairing order parameter oscillates as 
a function of position in the F side [s^]- The supercon- 
ductivity induced suppression of local density of states 
(LDOS) in the subgap energy range becomes enhanced 
where the order parameter changes its sign [3l| . Simi- 
larly, the internal tt state is expected to exhibit the zero 



2 



bias conductance enhancement in the subgap region [33| ■ 
We win show that it is indeed the case. This ehminates 
the major obstacle to estabhsh the sir state as the pair- 
ing symmetry of the FeAs superconductors. We solve the 
Bogoliubov-de Gennes (BdG) equation self-consistently 
to calculate LDOS of the tt state superconductor /normal 
metal (S/N) bilayers. It will be demonstrated that the 
observed zero bias conductance peaks which seemed diffi- 
cult to understand with a fully gapped pairing state may 
be consistently understood within the sir pairing state. 

Model - The FeAs superconductors have the discon- 
nected electron and hole Fermi surfaces. A minimal 
model has to include two bands; a hole band around the 
F = (0,0) point and an electron Fermi surface around 
the M = (7r,7r) point [s^l. We write the BdG equation 
for S/N bilayers with two band superconductivity as 
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Here, the subscript n ~ \ and 2 refer to the hole and 
electron Fermi surfaces around the F and M points in 
the momentum space, respectively, k is the intra-layer 
crystal momentum in the z-x plane, y is perpendicular 
to the interface between S and N. ^5 and are thickness 
of S and N layers, respectively, in the unit of the distance 
between neighboring single layers, and yt — ys ~^ yn- 1, 
ci, and C2 are the electron operators for the N, the hole 
band, and the electron band of S, respectively. The self- 
consistency relation for the gap function A„/c (y) is given 
by 

n' ,k' 

where V{n, k; n' , k') = Vnn' is the pairing interaction. 

Let us first consider bulk s-wave two band BCS super- 
conductivity. Putting VijNj — Xij, where iVi and N2 are 
the DOS per spin at the Fermi level for the band 1 and 
2, respectively, Eqs. ^ and ^ are reduced to [iol. Ib^I 



A2 = A2lF(Ai)Ai + A22i^(A2)A2, 

where F{A) is given by 

F(A)-^"°de ;^tanhQ/3i?^ , = V^M^. (7) 

At T = Tc, F = \n{l.UujD/Tc). The is the high- 
est temperature where the larger eigenvalue of Eq. Q 
becomes 1. 

To see physics through more clearly, consider the sim- 
plest case of All = A22 = A and A12 = A21 = A'. The 
single band BCS expression for the critical temperature 
Tc — I.IAujdc^^/^ is now replaced by one of the two 
expressions: 



rc = 1.14c^De-i/(^+^'), 



(8) 
(9) 



Ai = AiiF(Ai)Ai -f- Ai2i^(A2)A2 



(6) 



For A' < 0, Eq. ([9]) is the appropriate expression, and the 
pairing order parameter Ai and A2 on the two bands 
acquire the tt phase shift, which is the case considered 
here. The negative pairing interaction in one band BCS 
theory does not permit superconductivity. For the two 
band case, however, the negative interaction is turned to 
induce pairing by generating the sign reversal between Ai 
and A2 as can easily be seen from Eq. ([6]). The physical 
nature of the negative interaction parameterized in terms 
of A' we do not specify here, although it is most likely due 
to the antiferromagnetic fluctuations with a peak around 
the momentum transfer Q = (7r,7r) [13, 13, [Hi Hi]. 

The simple result of | Ai | = | A2 1 is an accidental conse- 
quence of the simple parameterization of An = A22 and 
A12 — A21. For more realistic parameterizations of the 
pairing interaction, the magnitudes of the two gaps are 
different and additional peaks show up in LDOS as shown 
in Fig. 2 below. One of many interesting consequences of 
the STT state is that the magnitude of one gap in general is 
larger than the BCS value while the other gap is smaller, 
that is, 2|Ai|/Tc > 3.52 and 2|A2|/7; < 3.52. An- 
other, perhaps more interesting feature of the stt pairing 
state is the appearance of the zero bias peak which may 
be probed experimentally by the tunneling spectroscopy. 
Szabo et al. recently performed directional point contact 
Andreev reflection spectroscopy on (Bao.55Ko.45)Fe2As2 
and found that some of the ab plane spectra reveal the 
zero-bias conductance peak consistent with the present 
work 37 1 . 

LDOS of S/N bilayers - Now, consider the S/N bi- 
layers of the STT pairing state described by Eq. ^ of the 
thickness yt — ys + Vn- The Hamiltonian is written as 
an M X M matrix, where M — Ays + '^yw, on the basis 
of which is taken as 

= (4^(1), Cl,-fei(l), 4^(1), C2,-fci(l), (10) 

cLt(2): ci,-fcx(2), • • • , al^{yt), akiiyt)) ■ 
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We first took the simple parameterization of A = 0, 
A' = —0.8, to ~ ti = t2 = 0.25 in the unit of the Fermi 
energy Ep. With the parameterization, we diagonahzed 
the M X M matrix and calculated the gap function us- 
ing Eq. (O. This procedure was repeated until the self- 
consistency was reached. In Fig. 1(a) we show the zero 
temperature 3 dimensional perspective plot of the LDOS 
of an S/N bilayer of 20 S layers and 20 N layers as a 
function of energy in the unit of the bulk pairing ampli- 
tude, V/Aq- The ratio of the gap to the Fermi energy 
in bulk, Aq/Ef, is 0.053. The coherence length in the 
unit of inter-layer distance is then ^5^3 — 4 and the 
thickness of S layer is large enough to monitor the evolu- 
tion of the proximity effects as one moves away from the 
S/N interface. Notice the LDOS enhancement in the sub- 
gap energy region near the interface. The origin of this 
zero bias enhancement is the sign change of the order pa- 
rameter. This point becomes clearer when we compare 
the present results with a two gap superconductor of the 
same sign. 




FIG. 1: (a) A 3D plot of LDOS of an stt state S/N bilayer 
with Ai — — A2 as a function of V/ Aq and the layer index. 
Notice the enhancement of the DOS in the subgap region near 
the interface due to the sign reversal of the order parameter, 
(b) The same as (a) but with Ai = A2. The LDOS behaves 
exactly as the well established conventional S/N bilayers. 

For comparison we repeated the same calculations as 
(a) with all the parameters remain unchanged except that 



the interband pairing interaction is turned to a positive 
A' = 0.8, corresponding to the case of Eq. ([8]). In this 
case, Ai — A2, and the LDOS should show the usual S/N 
behavior as shown in Fig. 1(b). The expected proximity 
effects are seen in the S and N regions. Compare this with 
the much short ranged proximity effects in the N layers 
of Fig. (a). The proximity effects of Ai and A2 for the sir 
pairing state cancel each other almost exactly other than 
the zero bias enhancement around the interface of S/N 
bilayers because they have the opposite sign and equal 
amplitude. 

The subgap enhancement is a robust feature of an stt 
state. It is a manifestation of the phase shift of tt between 
the two condensates which is insensitive to parameters. 
To demonstrate this we show in Fig. 2 the results of more 
realistic parameterization. We took An = 0.1, A22 = 0.2, 
A12 = -1.0, A21 = -0.5, to = 0.25, and ti = t2 = 0.2. 
As discussed above, other than another gap feature shows 
up inside the larger gap, the zero bias enhancement can 
clearly be seen. For the S/F tt state, the subgap en- 
hancement was observed for appropriate thickness of F 
where the order parameter changes its sign, referred to 
also as DOS reversal [sij. For the d-wave pairing case, 
the sign change of the order parameter occurs along the 
(110) surface, and the zero bias conductance peak was 
also pre dicted and observed in the cuprate superconduc- 
tors |38l . [sol , lioj . For stt state the peak appears unless 
one of ti and t2 is negligible such that both condensates 
are probed. 

Summary and outlook - We presented the local den- 
sity of states of an S/N bilayer based on a generic two 
band superconductivity model with the internal tt phase. 
The sign change of the pairing order parameter induced 
the LDOS enhancement in the subgap region near the in- 
terface as was observed by various tunneling spectroscopy 
on the FeAs superconductors. Some of the experiments, 
however, exhibit much sharper conductance peaks. This 
zero bias peak may be understood by more realistic ex- 
tensions of the present work. Dictated by the Hamilto- 
nian of Eq. ([T|) , an electron of the intra-layer momentum 
k tunneled from the N side may encounter the hole and 
electron band pairing gaps with the amplitudes ti and t2, 
respectively. More realistically, however, the specularly 
reflected electron and Andreev reflected hole may pick 
up the coherent phase difference of tt if their intra-layer 
momenta fall on the electron and hole Fermi surfaces, 
respectively. This process, which has been included on 
an average manner in the present work, can produce a 
much sharper conductance peak as observed in the pnic- 
tides and along the (110) surface of the cuprates. 

Other realistic considerations will be to allow the hop- 
ping amplitudes of the electron and hole bands with the 
N layers, ti and t2, to be different, or to consider more 
realistic pairing interactions. When these extensions are 
included, the LDOS will show more diverse and intrigu- 
ing behavior. Also interesting will be the spectroscopy of 
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FIG. 2: The LDOS of an tt state S/N bilayer with [Ail / 
I A2 1 . The blue and red curves are LDOS in the S and N layers, 
respectively. Each curve represents the LDOS of each layer 
and is shifted upward for clarity. Notice the enhancement of 
the DOS in the subgap region near the interface due to the 
sign reversal of the order parameter. 

other kinds of tt state multilayers. For instance, the S/F 
bilayers with the stt state will exhibit intriguing interplay 
between the internal and external vr phases. The works 
on these topics are in progress. 

Acknowledgement - This work was supported by Ko- 
rea Science & Engineering Foundation (KOSEF) through 
the Basic Research Program Grant No. ROl-2006-000- 
11248-0 (HYC) and by Korea Research Foundation 
(KRF) through Grant No. KRF-2005-070-C00044 (HYC, 
YB) and KRF-2007-521-C00081 (YB). 



[1] Y. Kamihara et al, J. Am. Chem. Soc. 128, 10012 
(2006). 

[2] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, 

J. Am. Chem. Soc. 130, 3296 (2008). 
[3] Phys. Today 61, Issue 5, 11 (2008). 
[4] G. F. Chen, Z. Li, D. Wu, G. Li, W. Z. Hu, J. Dong, P. 

Zheng, J. L. Luo, and N. L. Wang, Phys. Rev. Lett. 100, 

247002 (2008). 

[5] X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen, and D. 

F. Fang, Nature 453, 761 (2008). 
[6] C. de la Cruz, Q. Huang, J. W. Lynn, J. Li, W. Ratcliff 

II, J. L. Zarestky, H. A. Mook, G. F. Chen, J. L. Luo, N. 

L. Wang, P. Dai, Nature 453, 899 (2008)^ 

[7] M. Rotter, M. Tegel, and D. Johrendt, 'arXiv:0805.4630l 

(2008). 

[8] H. Ding, P. Richard, K. Nakayama, T. Sugawara, T. 



Arakane, Y. Sekiba, A. Takayama, S. Souma, T. Sato, 

T. Takahashi, Z. Wang, X. Dai, Z. Fang, G. F. Chen, 

J. L. Luo, and N. L. Wang, Europhys. Lett. 83, 47001 

(2008). 

[9] L. Zhao et aL. larXiv:0807.0398l 
[10] T. Kondo, A. F. Santander-Syro, O. Copie, C. Liu, M. 

E. Tillman, E. D. Mun, J. Schmalian, S. L. Bud'ko, 

M. A. Tanatar, P. C. Canfield, and A. Kaminski, 

larXiv:08 07.0815 
[11] C. Martin, R. T. Gordon, M. A. Tanatar, M. D. Van- 

nette, M. E. Tillman, E. D. Mun, P. C. Canfield, V. G. 

Kogan, G. D. Samolyuk, J. Schmalian, and R. Prozorov, 

arXiv:0807.0876 
[12] G. Li, W. Z. Hu, J. Dong, Z. Li, P. Zheng, G . F. Chen, 

J. L. Luo, and N. L. Wang. larXiv:0807. 10941 
[13] K. Matano, Z.A. Ren, X.L. Dong, L.L. Sun, Z.X. Zhao, 

and G.-q. Zheng, arXiv:0806.0249 
[14] H.-J. Grafe, D. Paar, G. Lang, N. J. Curro, G. Behr, J. 

Werner, J. Hamann-Borrero, C. Hess, N. Leps, R. Klin- 

geler, and B. Buechner, arXiv:0805.2595 
[15] Y. Nakai, K. Ishida, Y. Kamihara, M. Hirano, and H. 

Hosono, larXiv:0804.4765] 
[16] H. Mukuda, N. Terasaki, H. Kinouchi, M. Yashima, 

Y. Kitaoka, S. Suzuki, S. Miyasaka, S. Tajima, K. 

Miyazawa, P.M. Shirage, H. Kito, H. Eisaki, and A. lyo, 

arXiv:0806.3238 
[17] A. D. Christianson, E. A. Goremychkin, R. Osborn, 

S. Rosenkranz, M. D. Lumsden, C. D. Malliakas, 1. S. 

Todorov, H. Glaus, D. Y. Chung, M. G. Kanatzidis, R. 

I. Bewley, and T. Guidi, arXiv:0807.3932 
[18] T. Y. Chen, Z. Tesanovic, R. H. Liu, X. H. Chen, and C. 

L. Chien, Nature 453, 1224 (2008). 
[19] G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Phys. 

Rev. B 25, 4515 (1982). 
[20] L. Shan, Y. Wang, X. Zhu, G. Mu, L. Fang, and H.-H. 

Wen, arXiv:0803.2405 
[21] Y. Wang, L. Shan, L. Fang, P. Cheng, C. Ren, and H.-H. 

Wen, arXiv:0806.1986 
[22] O. MiUo, I. Asulin, O. Yuh, I. Felner, Z.-A. Ren, X.-L. 

Shen, G.-C. Che, and Z.-X. Zhao, ar Xiv:0807 .0359 
[23] P. Samuely, P. Szabo, Z. Pribulova, M. E. Tillman, S. 

Bud'ko, and P. C. Canfield, arXiv:0806.1672 
[24] I. I. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, 

Phys. Rev. Lett. 101, 057003 (2008). 
[25] D. Parker, O. V. Dolg ov, M. M. Korsh unov, A. A. Gol- 

ubov, and I. I. Mazin, larXiv:0807.3729l 
[26] A. V. Chubukov, D. Efremov, and I. Eremin, 

larXiv:0807. 3735 
[27] K. Kuroki, S. Onari, R. Arita, H. Usui, Y. Tanaka, H. 

Kontani, and H. Aoki, arXiv:0803.3325 
[28] A. G. Aronov and E. B. Sonin, Zh. Eksp. Teor. Fiz. 63 

1059 (1972) [Sov. Phys. JETP 36, 556 (1973)]. 
[29] M. J. Rice, H. Y. Choi, and Y. R. Wang, Phys. Rev. B 

44, 10414 (1991). 
[30] E. A. Demler, G. B. Arnold, and M. R. Beasley, Phys. 

Rev. B 55, 15174 (1997). 
[31] T. Kontos, M. Aprili, J. Lesueur, and X. Grison, Phys. 

Rev. Lett. 86, 304 (2001). 
[32] The zero bias peak due to the order parameter sign 

change has a deep root founded on the Atiyah-Singer 

index theorem. 

[33] S. Raghu, X.-L. Qi, C.-X. Liu, D. Scalapino, and S.-C. 

Zhang, Phys. Rev. B 77, 220503(R) (2008). 
[34] H. Suhl, B. T. Matthias, and L. R. Walker, Phys. Rev. 



5 



Lett. 3, 552 (1959). 
[35] Y. Bang and H.-Y. Choi, larXiv:0807.3912l 
[36] The sign reversal between the condensates means that 
the interband pairing must be repulsive, that is, 
Ai2, A21 < 0. This observation eliminates the phonon 
and leaves the Coulomb repulsion of both charge and 
spin channels for the pairing interaction. Together with 
this, the observations that the two bands are separated 
by the Q = (ti", tt) and the spin susceptibility has a peak 
also around Q strongly hint that the antiferromagnetic 



fluctuations act as the pairing glue. 
[37] P. Szabo, Z. Pribulova, G. Pristas, S. L. Budko, P. C. 

Canfleld, and P. Samuely, |arXiv:0809.1566 
[38] C. R. Hu, Phys. Rev. Lett. 72, 1526 (1994). 
[39] J. Y. T. Wei, N.-C. Yeh, D. F. Garrigus, and M. Strasik, 

Phys. Rev. Lett. 81, 2542 (1998). 
[40] H.-Y. Choi, Y. Bang, and D. K. CampbeU, Phys. Rev. B 

61, 9748 (2000). 



